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Abstract. We study the reionization histories where ionizing UV photons are 
emitted from decaying particles, in addition to usual contributions from stars and 
quasars, taking account of the fact that the universe is not fully ionized until z ^ 6 as 
observed by Sloan Digital Sky Survey. Likelihood analysis of the three- year data from 
the WMAP (Wilkinson Microwave Anisotropy Probe) severely constrains the decaying 
particle scenario. In particular, the decaying particle with relatively short lifetime is not 
favored by the polarization data. 
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1. Introduction 

Ionization history after recombination was not known for many years, until recent 
development of observations. The Sloan Digital Sky Survey (SDSS) obtained the first 
evidence of the reionization epoch [1]. The Gunn- Peterson test for the observed spectra 
of quasars revealed that there existed neutral hydrogen at the redshift 2; ~ 6, which 
implies that the completion of the reionization of the universe was seen at 2; ~ 6, 
although the amount of the neutral hydrogen aX z > 6 may not be very large to explain 
the observed Gunn-Peterson trough. 

On the other hand, the Wilkinson Microwave Anisotropy Probe (WMAP) 
observation found the earlier reionization than ^ ~ 6 in terms of the large optical 
depth [21 [3]. The first-year data favored very large optical depth Top = 0.17 ± 0.04 [2]. 
The conventional reionization by UV photons emitted from usual early formed stars and 
quasars, could only reach to the optical depth of ~ 0.05 [H [5] if it is consistent with 
the SDSS data. Of course, using the nonstandard initial mass function and unrealistic 
escaping UV fraction from galaxies, one could achieve a large optical depth [6]. However, 
the UV fluxes from stars and quasars increase with time and full reionization takes place 
before z ~ 6. Thus, one has to consider a more complex ionization history which may 
be realized, for example, by the use of different photon emission processes between 
Population II and III stars [7j. 

An alternative line is to consider another source of UV photons. It could be decaying 
particles [8], [9l [10]. In Refs. [8], [9], we considered decaying particles in addition to 
conventional stars and quasars and found that the decaying particle scenario is consistent 
with the first year WMAP data for large parameter space. 

Recently, WMAP team released three year data, in which the value of the optical 
depth is lowered: Top = 0.09 ± 0.03 [3]. It changed towards the conventional value of 
~ 0.05, but some extra UV photons are still necessary. It would be possible to have such 
a small amount of extra UV photons before z ~ 6 using less unrealistic assumptions for 
stars and quasars. In this article, however, we stick to the decaying particle scenario to 
assess its capability. 

The structure of the paper is as follows. In the next section, we explain how 
we follow the evolution of the ionization history including the decaying particles. In 
Sect. Ill, we evaluate the for each model parameters, and construct likelihood 
contours for the parameters of decaying particle scenario. At the same time, the 
detail of the power spectra of temperature and polarization anisotropics in the cosmic 
microwave background (CMB) radiation are shown, and we discuss how it discriminates 
the preferred model parameters in the same section. The final section is devoted to our 
conclusions. 
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2. Reionization history 

We consider that ionizing UV photons are emitted from decaying particles as well as 
from usual stars (and quasars). The latter are indeed responsible for full ionization at 
2; ~ 6. On the other hand, UV photons from particle decays keep a nonzero but small 
amount of the ionization fraction before 2; ~ 6. 

We follow the thermal history from z > 10^ including the recombination epoch, 
calculating the ionization fraction of hydrogen and helium, and the electron temperature, 
on the basis of the hierarchical clustering scheme of the cold dark matter scenario used 
in Fukugita and Kawasaki fTT]. 

In addition, we include the contribution of UV photons from decaying particles 
[H [9]. We assume that the particle emits two photons with monochromatic energy 
of half the mass of that particle, i.e., = The number density of (p particle is 

given by 

n^ = r^^(0)(l + 2)V*/^^ (1) 

where is the lifetime of the particle. Hydrogen atoms can be ionized if the emitted 
photons have energy with E^ > 13.6 eV. Then the source term for the decaying particle 
can be written as 

The amount of emitted photons is determined once the mass, m^, lifetime, r^, and 
abundance, fi^, of the particle are fixed. We calculate the ionization histories for 
E^ = 15 — 1000 eV and = 10^^ — 10^^ sec, adjusting the abundance, fl^, to get 
the optical depth in the range of Top = 0.075 — 0.2, which is defined by 

dzarl^] K-ni-)], (3) 



where ax is the Thomson coss section and n^^"^^ is the electron number density 
for standard recombination. We subtract this term to estimate only the effect of 
reionization. 



3. Cosmic microwave background radiation 



As well known, the temperature and polarization anisotropics of the CMB radiation 
could inform us of what happened during the dark age after recombination. We now 
look for whether ionization histories in the decaying particle scenario are consistent with 
CMB observation in the three-year WMAP data. 

We obtain the power spectra by the code modified from CMBFAST [12] so as to 
equip us with the capability of using the evolution obtained for the ionization fraction 
and matter temperature, and evaluate using the code provided by WMAP [3j. We 
seek for the range E^ = 15 — 1000 eV, = lO^'' — 10^^ sec (with 9 bins), and 
Top = 0.075 — 0.2 (with 6 bins). For fixed mass, lifetime, and optical depth, is 
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calculated adjusting fi^, Vtm, Hq, Us, and the amplitude of the spectrum. We construct 
the contour of likehhood with grid-based analysis, using Ax^ = X^^X^Ibest, where x^lbest 
is the six-parameter minimum for ACDM with WMAP-only. 

Before going into the analyses with hybrid UV source of stars and decaying particles, 
let us ask a question: Is it possible that the decaying particles alone could provide an 
ionization history consistent with SDSS and three-year WMAP? It is indeed possible to 
have an ionization history such that the ionization fraction rises rapidly at z ~ 6, which 
is consistent with SDSS observation. For example, this is realized with = 9.5 x 10^^ 
sec, Q(f, = 1.2 X 10^^, and = 80 eV. However, the is about Aa larger than that of 
the best-fit ACDM model. Of course, the optical depth has a very large value, such as 
0.23. Thus, we can exclude this possibility. 

We next consider another limiting case, that the ionization fraction becomes unity 
at 2; ~ 6 due to UV photons emitted from stars and quasars only, and there is no 
contribution from decaying particles, which leads to the optical depth Top ~ 0.05. In 
this case, the is just outside the la region. One may hope that the will be reduced 
once an extra contribution of UV photons from decaying particles is included. As we 
shall see below, however, this is not the case. 

We show the likehhood contours for = 15, 100, and 1000 eV, in the top-left, 
top- right, and bottom panels, respectively, in Fig. [H They look very similar to each 
other. Contrary to our expectation, there is no Icr region. We can see that there is a 
weak dependence on the lifetime, and that the larger optical depth is not favored. 

It is wise to compare these with the contours for which only the temperature power 
spectrum is used for estimating the likelihood. We can then see that the does not 
depend on the lifetime at all, as can be seen in Fig. [2] for E^ = 100 eV. Moreover, 
there is a la region: Top < 0.13(0.17) is allowed for the 68% (90%) confidence level. 
This is easy to understand, because the temperature power spectrum is dependent on 
only the value of the optical depth, and does not 'see' the actual ionization history. On 
the other hand, the polarization anisotropy in the CMB radiation discriminates among 
the different ionization histories. This is why we can see the dependence on the lifetime 
in Fig. [H and the longer lifetime seems to be preferable. 

Notice that lifetime longer than 10^^ sec for E^ = 1 keV is excluded. This is because 
the flux of the photons exceeds the observed diffuse X-ray background p], although the 
ionization history respects the Gunn-Peterson trough observed by SDSS. 

We can thus obtain the allowed region for the decaying particle scenario in the 
parameter space (r^,^^). We show each line with a 99% confidence level in Fig. [31 
For Ej < 100 eV, longer lifetime is also allowed, until the abundance reaches the whole 
amount of the dark matter, since it is the ratio of u^/t^ that determines the amount 
of the emitted photons from the decaying particles for lifetime longer than the age of 
the universe. Among the decaying particles, saxion, the scalar partner of the axion, 
in the gauge-mediated supersymmetry breaking scenario, may be the best candidate. 
Its lifetime is about 10^^ sec for the saxion mass iris ~ 100 eV and the axion decay 
constant Fpq ~ 10^^ GeV. Since the saxion decays into photons through a one-loop 
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Figure 1. Likelihood contours of the decaying particle scenario with = 15 eV 
(top-left), 100 eV (top-right), and 1 keV (bottom). We show 90, 95, and 99% CL 
lines. Lifetime longer than ~ 10^^ sec (shaded region) is excluded by diffuse X-ray 
background observation in the bottom panel. 
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Figure 2. Likelihood contours of the decaying particle scenario with ~ 100 eV, 
using only temperature anisotropy data. We show 68 and 90% CL lines. 



diagram with branching ratio of order ~ 10 — 10 ^, it matches to the allowed region 
of the parameter space very well. Notice that these constraints apply to most kinds 
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of decaying particles which emit photons with branching ratio provided that one 
regards the effective abundance as B^VL^. 
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Figure 3. Allowed region of the parameter space of the decaying particle scenario. 
The lines are drawn for the 99% CL contours in Fig.[Tl which are thus the upper limits 
for the amount of the decaying particles, Sl^, for fixed lifetime r^. 

Now let us move on to the detailed power spectra of temperature and 
polarization anisotropics in CMB. Temperature (TT), temperature-polarization (TE), 
and polarization (EE) power spectra are shown for optical depth in the range of 
Top = 0.075 — 0.2 with fixed lifetime = 10^^ sec in Fig. HI For the limiting case 
of the decaying particle scenario, we also plot the spectra without any contributions 
from decaying particles {jop ~ 0.05). In addition, one-step instantaneous reionization 
case ('?~op 

~ 0.09) is shown. 

Since it is not easy to obtain some information only from the figures of the spectra, 
we give the breakdown of the in these cases in upper half of Tabled! Here, we tabulate 
the difference of from the minimum case, Ax^, defined above. Larger optical depth 
is disfavored by all the components, which may be seen in all the spectra in Fig. HI The 
case with Top ~ 0.01 seems to be slightly favored by the polarization data, which almost 
mimic the instantaneous reionization case in TE spectra. However, the significance is 
not so great, and TT data do not prefer this case. Anyway, the total slightly favours 
smaller optical depth. 

On the other hand, one can look for the dependence on the lifetime with fixed 
optical depth in Fig. [5] and the lower half of Table [H In this case, TT data do not 
favour any lifetime, since the TT power spectra is almost only sensitive to the value 
of the optical depth, not to the history of actual ionization. Short lifetime is severely 
disfavored since it has large x^ from the polarization data. This is not so apparent from 
the figures for either TE or EE power spectra, but it seems that larger power is needed 
in the low £ regions in order to fit the polarization data much better. 
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Figure 4. TT (top-left), TE (top-right), and EE (bottom) power spectra of the 
CMB for Tff, = 10^^ sec with various optical depths. We also show the spectra for the 
limiting case of no decaying particle contribution, and the (one-step) instantaneous 
reionization. The three year WMAP data are plotted as well. 



4. conclusions 

We have considered the ionizing UV photons emitted from decaying particles, in addition 
to the usual contributions of stars and quasars. The stars and quasars are responsible 
for the quick rise of the inonization fraction at z ~ 6 observed by SDSS as the Gunn- 
Peterson trough. On the other hand, the decaying particles cause partial ionization at 
z>6. 

The decaying particle scenario could explain the large optical depth, and was 
consistent with the first-year WMAP data. However, the three-year WMAP data do 
not favour the decaying particle scenario so much. Notice that our results apply to most 
kinds of decaying particles which emit photons with branching ratio provided that 
one regards the effective abundance as B^Q^. 

The ionization by the decaying partcles with Top < 0.17 does indeed seem to be 
consistent within 2a with TT power spectrum. The peculiar aspect of the three-year 
WMAP data is in the polarization; in particular, the polarization data do not allow 
the earlier reionization which is a distinct feature of the ionization history from the UV 
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Table 1. Ax^ for various optical depth for T0 = 10^^ sec (upper), and for various 
lifetimes for Top w 0.1 (lower). In the upper half, the last line represents for the case 
without contributions from decaying particles. TT, low TT, TE, and low pol. stand 
for the TT power spectrum, pixel based TT for low £, TE power spectrum, and pixel 
based polarization for low £, respectively. 
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1.16 


0.01 


1.42 


2.65 


5.25 


1014 


0.100 


1.21 


-0.23 


4.04 


2.79 


7.82 



emission from particle decays, and hence the longer lifetime is favoured. Since more 
power at low i is favoured by the polarization anisotropy, one may ambitiously claim 
that a better fit may be obtained for the decaying particle scenario with increasing decay 
rate. 
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(2003); A. Sokasian, N. Yoshida, T. Abel, L. Hernquist and V. Springel, Mon. Not. Roy. Astron. 



Soc. 350, 47 (2004); A. J. Benson, N. Sugiyama, A. Nusser and C. G. Lacey, |astro-ph/0512364 
[8] S. Kasuya, M. Kawasaki and N. Sugiyama, Phys. Rev. D 69, 023512 (2004). 
[9] S. Kasuya and M. Kawasaki, Phys. Rev. D 70, 103519 (2004). 
[10] S. H. Hansen and Z. Haiman, Astrophys. J. 600, 26 (2004); X. L. Chen and M. Kamionkowski, 

Phys. Rev. D 70, 043502 (2004); P. P. Avehno and D. Barbosa, Phys. Rev. D 70, 067302 (2004); 

P. L. Biermann and A. Kusenko, Phys. Rev. Lett. 96, 091301 (2006); M. Mapelh, A. Ferrara and 

E. Pierpaoh, Mon. Not. Roy. Astron. Soc. 369, 1719 (2006); L. Zhang, X. L. Chen, Y. A. Lei 



and Z. G. Si, astro-ph/0603425; Y. A. Shchekinov and E. O. Vasiliev, astro-ph/0604231 



[11] M. Fukugita and M. Kawasaki, Mon. Not. R. Astron. Soc. 269, 563 (1994). 
[12] U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 (1996). 



